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SUMMARY 


We  develop  a  simple  analytical  model  to  estimate  the  thickness  of 
a  smoke  layer  formed  by  the  plume  of  a  large  area  fire  and  to  account 
for  crosswinds.  We  take  advantage  of  the  dominant  flow  features  in 
the  upper  part  of  the  rising  plume  and  in  the  smoke  layer  far  from  the 
plume  to  model  the  transition  of  a  vertical  flow  to  a  horizontally 
spread  thin  smoke  (or  dust)  layer.  A  one-dimensional  approximation  is 
used  to  model  the  buoyancy  decay  and  plume  overshoot.  That  solution 
is  matched  and  combined  with  a  solution  based  on  the  work  required  to 
establish  our  asymptotic  horizontal  smoke  layer.  The  solution  es¬ 
timates  the  thickness  of  the  spreading  smoke  layer  in  a  quiescent 
atmosphere  or  in  an  ambient  wind  field. 
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PREFACE 


This  report  continues  Pac i f ic -Research  Corporation's  study  of 
nuclear  weapon  fire  effects.  In  this  volume  we  develop  a  simple  model 
to  predict  the  stabilization  altitude  and  thickness  of  smoke  layers 
formed  by  large  fires.  Such  smoke  layers  spread  over  large  distances 
and  can  obscure  large  areas.  The  cloud  altitude  and  thickness  are 
inputs  for  meso-  and  global-scale  environment  calculations. 

This  research  was  supported  by  the  Defense  Nuclear  Agency  under 
contract  DNA  001-88-C-0119  and  was  monitored  by  Dr.  David  Auton,  RARP. 
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SECTION  1 
INTRODUCTION 


The  spread  of  smoke  over  regional  scales  depends  on  the  initial 
smoke  injection  and  the  temperature  and  wind  structure  of  the 
atmosphere.  Although  the  tran.sition  from  a  vertical  plume  flow  to  a 
biorizontal  cloud  is  three-dimensional,  separately,  the  spread  and  the 
injection  can  be  approximated  as  two-dimensional  flows.  Regional 
particulate  cloud  movements  are  often  in  thin  layers  with  only  slow 
changes  normal  to  the  principal  flow  direction  [Westphal,  et  al., 
1988'.  Similarly,  large  fire  plume  motions  are  principally  vertical 
and  two-dimensional  ( axisvmme tr ic )  [Small  and  Larson,  1985;  Small  and 
Heikes,  1988':.  A  similar  division  of  flows  also  occurs  when  the  dust 
pb-ime  following  a  rising  fireball  is  sheared  at  the  stabilization 
al t i tude . 

Plumes  formed  by  nuclear  weapon  fires  [Small,  1989;  can  have 
diameters  comparable  to  the  atmospheric  scale  height.  The  plume 
motion  i.s  thus  a  significant  perturbation  to  the  atmosphere.  Gravit;.’ 
waves  which  propagate  long  distances  are  formed  by  dynamical  processes 
in  the  fire  inflow  region  or  by  the  oscillation  of  the  plume  about  its 
neutral  buoyancy  level  in  the  upper  atmosphere  [Heikes  and  Small, 

1990a:  .  Large  fire  plumes  can  penetrate  high  into  the  atmosphere; 
their  rise  limited  mostly  oy  gradients  in  temperature  and  winds. 
Because  nuclear  weapon  fires  can  cover  several  hundred  square 
kilometers  and  plume  diameters  can  be  roughly  10  km  or  greater,  am¬ 
bient  winds  do  not  greatly  influence  the  plume  motion  over  most  of  its 
rise.  However,  near  the  top  of  the  plume,  the  buoyancy  is  reduced  and 
the  motion  slows,  upper  level  winds  bend  the  column  and  the  plume  is 
sheared  into  thin  clouds  that  are  advecteu  away  from  the  target  area. 
The  thickness  of  the  spreading  cloud  can  be  calculated  [Small  and 
Heikes,  1989;  Heikes  and  Small,  1990aj  in  quiescent  atmospheres.  Two- 
dimensional  or  axisymmetri'c  hydrocode  solutions  can  calculate  the 
pbume  rise,  overshoot,  and  settling  at  the  neutral  buoyancy  level. 

Such  solutions,  however,  do  not  account  for  crosswinds;  three- 
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dimensional  hydrocodes  are  required  and  only  a  few  solutions  have  been 
developed  [Cotton,  1985;  Bradley  et  al . ,  19^8].  In  those  solutions, 
the  resolution  is  fairly  coarse  and  the  cloud  thickness  far  from  the 
plume  is  only  estimated. 

An  alternative  to  three-dimensional  code  solutions  uses  modifica¬ 
tions  of  the  Morton  et  al.  [1956]  plume  similarity  solution  to  account 
for  cross-winds  [Briggs,  1975;  Manins ,  1985].  Such  solutions  are 
valid  for  asymptotically  slender  plumes  in  which  the  buoyancy  is 
diffused  by  entrainment  of  ambient  air.  For  nuclear  weapon  fire 
plumes--or  even  large  prescribed  fires  [Heikes  and  Small,  1990a] --the 
plume  is  too  thick  for  entrainment  at  the  edge  to  diffuse  the  buoyancy 
or  momentum.  Similarity  type  solutions  have  been  applied  to  plumes 
from  industrial  or  utility  smoke  stacks.  They  are  not  applicable  for 
nuclear  weapon  smoke  or  dust  plumes. 

In  this  report,  we  develop  a  simple  two-dimensional  model  to 
estimate  the  thickness  of  the  smoke  (or  dust)  cloud  far  from  the  plume 
and  to  account  for  crosswinds.  We  take  advantage  of  the  dominant  flow 
features  in  the  rising  plume  and  far  field  to  model  the  transition  of 
a  vertical  flow  to  a  horizontally  spreading  thin  smoke  (dust)  layer. 

A  one-dimensional  approximation  is  used  to  model  the  buoyancy  decay 
and  plume  overshoot.  That  solution  is  matched  and  combined  with  a 
solution  based  on  the  work  required  to  establish  a  asymptotic  horizon¬ 
tal  smoke  layer.  The  solution  while  approximate  estimates  the  thick¬ 
ness  of  clouds  that  merge  (for  multiple  plumes)  and  spread  regionally. 
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SECTION  2 


ANALYSIS 


In  the  upper  portion  of  the  plume  region,  the  combustion  products 
and  entrained  air  have  been  moving  upwards  for  a  time  of  order  z/W 
where  the  height  z  is  of  order  1000  m  and  the  average  vertical 
velocity  W  is  of  order  10  m/s.  Even  for  low  turbulence  levels,  the 
plume  is  well  mixed  and  in  approximate  thermal  equilibrium.  The 
density  deficit  (buoyancy),  pressure,  and  the  velocity  are  thus  well 
defined,  continuous  functions  of  position. 

The  full  flowfield  describing  the  formation  of  a  capping  cloud  is 
subtle,  involving  the  interaction  of  small  temperature,  density  and 
pressure  differences  within  and  external  to  the  plume.  Howe”er,  by 
breaking  the  flowfield  into  zones,  most  of  the  important  features  of 
the  upper  plume  can  be  described  by  simple  models.  These  include  the 
bending  of  the  plume  axis  by  crosswinds,  establishing  the  maximum 
height  obtained  by  the  combustion  products,  and  determining  the  shape 
and  thickness  of  the  outward  spreading  smoke  layer  as  a  function  of 
distance  from  the  plume  center.  The  altitude  and  thickness  of  the 
initial  layer  is  the  initial  condition  for  calculating  the  mesoscale 
spread  of  smoke.  The  cloud's  position  relative  to  the  tropopause 
influences  the  rate  of  scavenging.  At  low  altitudes  precipitation 
enhances  particle  removal  while  in  the  stratosphere  removal  is  slow. 

The  fire  induced  flowfield  can  be  described  by  six  zones,  each  of 
which  is  dominated  by  different  physical  mechanisms.  Simple  equation 
sets  can  be  separately  formulated  for  each  zone,  and  the  results 
patched.  The  zones  are  (see  Fig.  1): 


1.  Inflow  layer.  Approximated  by  a  sink-like  flowfield  [Weihs 
and  Small,  1986].  In  this  zone  air  moves  at  low  (incompres¬ 
sible)  speeds  close  to  ground  level.  The  inflow  (fire  winds) 
provides  oxygen  for  the  fire;  it  is  driven  by  t}ie  pressure 
deficit  in  the  burning  zone.  Here  r  >  Rf,  z  <  rig. 

2.  The  flame  and  combustion  zone.  Defined  as  a  cylindrical 
volume  of  radius  Rf  and  height  Hf  in  which  exothermic  reac¬ 
tions  take  place;  combustibles  and  air  are  turned  into  hot 
buoyant  plume  gases. 


3 


Figure  1.  Schematic  representation  identifying  principal  flov/  regions 
of  axi symmetric  plume  driven  by  a  large  area  fire. 


3.  Turning  region.  Here  the  radial  inflow  reduces  the  diameter 
of  the  upward  moving  plume.  Initially,  the  contraction  angle 
is  quite  acute.  The  slope  is  reduced  at  higher  elevations  as 
the  plume  contracts,  accelerating  the  vertical  motion,  while 
inflow  is  reduced.  The  vertical  bounds  of  zone  3  are  Hf  <  z 
<  Hq.  Buoyance  is  initially  high,  and  heat  loss  by  radiation 
is  a  significant  factor. 

4.  Central  column.  This  is  a  major  part  of  the  convection 
column;  the  flow  is  essentially  parallel,  moves  upwards  and 
can  be  approximated  as  one-dimensional.  Temperature  dif¬ 
ferences  and  buoyant  forces  are  small.  Shear  at  the  column 
perimeter  affects  only  a  relatively  thin  boundary  layer.  At 
the  upper  end  of  zone  4,  the  plume  density  is  equal  to  that 
of  the  atmosphere,  and  crosswinds  start  shearing  the  plume. 

5.  Capping  cloud.  In  this  region  the  plume  density  is  essen¬ 
tially  constant  and  greater  than  the  surrounding  air.  The 
upward  momentum  causes  the  plume  to  rise  further  and  finally 
fall  back,  while  crosswinds  advect  it  downstream.  It  then 
spreads  horizontally.  The  flow  in  this  zone  is  intrinsically 
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Outflow 


three-dimensional,  ending  at  a  distance  where  the  vertical 
component  is  much  smaller  than  the  horizontal  one.  This 
defines  a  critical  radius  (R^r) • 

6.  Equilibrium  cloud  layer.  This  is  a  region  of  purely  horizon¬ 
tal  flow,  at  a  height  H*  .  This  zone  is  an  annular  disk  of 
radius  r  >  with  slowly  varying  width  (the  final  cloud 

layer  thickness)  distorted  by  ambient  crosswinds. 

2.1  CENTRAL  COLUMN  ANALYSIS. 

We  first  consider  the  flow  in  zone  4  without  ambient  winds.  This 
analysis  is  the  starting  point  for  including  crosswinds  and  estimating 
the  maximum  overshoot  height  of  zone  5). 

Zone  4  is  characterized  by  weak,  decreasing  buoyancy  and  negli¬ 
gible  entrainment.  Initial  conditions  at  Hq  include  the  vertical 
velocity  and  the  temperature  excess.  These  are  obtained  from  either 
analytic  solutions  [Small  and  Larson,  1985]  or  from  hydrocode  simula¬ 
tions  of  the  fire  zone  [Bacon  et  al.,  1987;  Penner  et  al . ,  1986;  Small 
and  Heikes,  1989;  Heikes  and  Small,  1990a, b] .  Maximum  speeds  are 
generally  0[10]  m/s  and  the  temperature  excess  is  small,  usually  less 
than  10  K.  The  largest  values  of  both  are  at  the  plume  axis  (this  may 
be  an  artifact  of  the  radial  symmetry,  although  in  large  experimental 
fires  [Radke,  1990],  the  plume  top  is  generally  near  the  fire 
center) , 

To  establish  the  dominant  flow  features,  we  make  several  assump¬ 
tions  that  simplify  the  formulation.  We  assume  steady  flow.  The 
analysis  is  thus  valid  after  the  initial  ignition-fire  plume  formation 
transients  are  over  and  almost  to  burnout.  The  initial  transients 
generally  last  5  to  20  minutes.  Typical  fires  last  1  to  4  h  [Small, 
1989]  so  that  over  80  percent  of  the  fire  duration  is  covered  by  the 
present  analysis. 

The  plume  is  assumed  mostly  air  with  properties  of  an  ideal  gas. 
Most  fires  burn  with  excess  air  and  the  combustion  products  generally 
comprise  less  than  several  percent  of  the  total,  even  if  perfect 
combustion  is  achieved  (stoichiometric  ratios  for  most  prevalent 
combustibles  are  around  15:1).  The  amount  of  water  added  by  combus¬ 
tion  products  is  usually  small  relative  to  entrained  ambient  humidity, 
except  in  particularly  arid  areas  [Small  and  Heikes,  1987]. 
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Heat  transfer  between  the  plume  and  surroundings  is  neglected. 
Zone  3  may  have  significant  heat  loss  due  to  radiation  [Murgai,  1962; 
Smith,  1967],  Above  that  region,  temperature  differences  are  small, 
and  adiabatic  or  pseudo-adiabatic  flow  is  assumed  in  zones  4,  5,  and 

6. 

Mass  transfer  between  the  plume  and  the  surroundings  is 
negligible.  This  is  equivalent  to  assuming  that  all  entrainment  takes 
place  in  the  lower  zones,  and  that  the  plume  mass  flux  in  zone  4 
changes  only  slightly.  Thus,  the  limiting  streamlines  (recalling  the 
assumption  of  steady  flow)  defined  at  the  lowest  plane  of  zone  4 
outlines  the  plume  trajectory. 

The  structure  of  the  external  atmosphere  is  not  strongly  in¬ 
fluenced  by  the  presence  of  the  plume,  and  the  boundary  conditions  for 
external  pressure,  density,  and  temperature  are  independent  of  time 
and  fire  parameters. 

The  altitude  H*  defining  the  equilibrium  height  of  the  cloud,  as 
well  as  the  upper  end  of  zone  4,  and  the  vertical  velocity  (W*)  at  the 
altitude  is  obtained  from  one  dimensional  flow  considerations.  The 
equation  of  continuity  now  reads  simply 

pwA  =  Const  ,  (1) 

where  A  is  fixed,  and  input  conditions  at  Hq  obtained  from  a  separate 
analysis  or  numerical  results  [Small  et  al.,  1988].  Equation  (1)  is 
rewritten  as 


* 

W 


(2) 


Equation  (2)  provides  the  value  W*  once  the  density  equation  is 
found.  The  density  changes  in  the  plume  are  adiabatic  so  that 


(3) 
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The  ambient  atmosphere  has  a  lapse  rate  dependent  on  time  of  day, 
season,  and  humidity.  This  requires  an  independent  set  of  functions 
to  describe  the  altitude  variation  of  pressure  and  density.  Addi¬ 
tional  information  is  required  to  find  the  equilibrium  density,  that 
is  the  density  where  both  Eq .  (3)  and  the  specific  meteorology  give 
the  same  values  for  both  pressure  and  density.  This  density  then 
defines  the  equilibrium  height  H*,  and  by  Eq .  (2)  the  velocity  W*. 

The  density  distribution  is  affected  by  water  vapor.  The  main  effect 
is  due  to  moisture  (ambient  humidity)  entrained  by  the  firewinds 
(water  produced  as  a  product  of  combustion  is  small  compared  to  the 
entrained  moisture).  For  the  following  calculations,  we  use  a  stand¬ 
ard  relative  humidity  model  [Manabe  and  Wetherall,  1967] 


RH(z) 


RH(0)  P(z) 
0.98  P(o) 


0.02 


The  humidity  at  the  bottom  of  zone  4  is  the  average  humidity  over  the 
height  range  0  <  z  <  1000m.  The  ambient  humidity  produces  saturation 
(except  in  very  arid  cases)  as  the  ambient  temperature  falls  and  the 
dew-point  is  crossed.  Most  of  the  condensation  occurs  within  zone  4. 
The  latent  heat  released  reduces  the  density  which  effectively  in¬ 
creases  K  in  Eq .  (3). 


2.2  ANALYSIS  OF  CAPPING  CLOUD  REGION. 

This  same  model  is  also  used  to  calculate  the  maximum  height 
obtained  by  the  plume  •  This  calculation  assumes  [Turner,  1973] 

that  density  is  constant  in  overshoot  regions  such  as  zone  5.  The 
maximum  height  is  achieved  by  the  central  streamline  so  that  a  one¬ 
dimensional  calculation  again  is  sufficient.  For  these  assumptions, 
the  momentum  equation  in  the  vertical  direction  reduces  to 


P 

pA^) 

a 


1 


(4) 
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where  p^{z)  is  the  ambient  density.  Integrating  Eq .  (4)  over  z  gives 
the  vertical  velocity  distribution  as  a  function  of  ambient  density 
and  the  initial  profile  from  zone  3. 

Results  of  these  computations  for  a  7  km  diameter  fire  and 
average  heat  addition  rate  of  100  KW/m^  are  shown  in  Figs.  2  and  3. 
Input  conditions  at  z  =  1000  m  are  AT  =  10  K  and  Wq  =  10  m/s  [Small 
and  Heikes,  1988].  The  standard  day  absolute  humidity  is  nominally 
7.5  g  H20/kg  dry  air.  For  this  fire,  the  condensation  height  was 
~2500  m,  the  equilibrium  smoke  layer  height  -11  km,  and  the  maximum 
plume  height  about  15  km.  The  comparison  of  predicted  cloud  altitudes 
in  Fig.  2  shows  that  the  present  approximation  satisfactorily  es¬ 
timates  the  maximum  altitudes  within  the  range  and  spread  of  current 
hydrocode  results. 

2.3  PLUME  RISE  IN  CROSSWIND. 

Next,  we  consider  the  plume  rise  in  zone  4  when  subjected  to 
ambient  crosswinds.  Annually  averaged  wind  velocities  as  a  function 
of  height  above  sea  level  are  given  in  Fig.  4.  The  "standard"  profile 
averages  wind  profiles  used  by  Penner  et  al.,  [1986]  and  Scoggins  and 
Vaughan  [1962].  The  latter  windspeeds  are  0[10  m/s]  (comparable  to 
the  upflow  speeds  in  zone  4) . 

To  enable  a  serai-analytic  solution  for  the  plume  deflection  by 
such  crosswinds,  several  simplifying  steps  are  made.  We  assume  that 
plume  horizontal  cross-sections  stay  round  and  that  the  rising  plume 
gases  do  not  mix  with  the  ambient  air,  i.e.,  the  plume  boundary  is  a 
streamtube.  Each  horizontal  cross-section  is  assumed  to  be  displaced 
relative  to  its  vertical  neighbors  under  the  influence  of  the  wind  at 
its  height.  The  resulting  bending  of  the  plume  can  thus  be  visualized 
as  a  series  of  disks  pushed  sideways  by  a  force  varying  with  height. 

There  are  three  possible  contributions  to  the  horizontal  force. 
These  are  the  aerodynamic  drag  force  between  wind  and  plume,  the 
centrifugal  force  resulting  from  plume  bending,  and  the  Coriolis  force 
which  results  from  the  angle  between  the  vertical  plume  and  the 
Earth's  rotation  vector  (complementary  angle  to  latitude).  The 
Coriolis  acceleration  is 
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a  =  2Qw  cos(fi  ,  (5) 

where  n  =  0.0073  rad/s  is  the  Earth's  angular  velocity,  w  is  the 
vertical  velocity  of  the  plume  and  4>  is  the  geographic  latitude.  We 
directly  obtain  the  horizontal  displacement  of  the  plume  centerline  as 
a  function  of  fluid  particle  "residence  time"  in  zone  4,  =  AH/w. 

This  is 

Ax  =  a^  r^/2  =  n  (AH/w)V2  ,  (6) 

where  AH  =  15000  m,  and  100  >  W  >  20  m/s,  so  that  <  100  s ,  and 

Ax  <  100  m  ,  (7) 
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i.e.,  much  smaller  than  the  other  effects,  as  we  shall  demonstrate. 

We  therefore  neglect  the  displacement  caused  by  the  Coriolis  effect  in 
the  following  analysis. 

The  other  two  forces  are  the  drag  on  the  plume  (the  pushing 
force)  and  the  centrifugal  force  resulting  from  the  plume  bending. 
These  are  equal  and  opposite,  in  the  direction  normal  to  the  curved 
plume.  The  drag  force  per  unit  plume  length  is 

dD  =  p  cos^d  •  b  C,  dz  ,  (8) 

aw  d 

where  p^  is  the  ambient  air  density,  is  the  windspeed,  9  is  the 
local  angle,  b  is  plume  radius,  and  is  a  drag  coefficient.  The 
centrifugal  force  obtained  in  the  curved  plume  is 

,  , 2  ,  2  d^ 

dz  =  pUjrb  —  dz,  (9) 

p  t  dz 

where  pp  is  the  plume  density,  Uj-  is  the  total  ("including  both  verti¬ 
cal  and  horizontal  components)  of  the  velocity  and  d9/dz  is  the  local 
change  of  plume  angle. 

From  Eqs .  (9)  and  (2),  we  obtain 


dF  =  m 


dU. 


dz 


cos 


dz 


(10) 


which  is  integrated  to  obtain  the  angular  deflection  of  the  rising 
plume  as  a  function  of  height.  Additional  integration  yields  the 
horizontal  motion  of  the  plume  centerline.  A  parametric  study  of 
pl^ne  centerline  drift  for  the  nominal  7  kill  fire  is  presented  in 
Fig.  5.  Three  values  are  chosen  for  the  drag  coefficient.  The 
highest  value  (0.7)  is  that  of  a  solid,  unmoving  cylinder  in  high 
Reynolds  number,  fully  turbulent  flow.  The  plume  boundary  is  not 
aolid,  so  that  the  crosswind  flow  will  cause  some  parallel  motion  in 
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Figure  5.  Plume  centerline  as  a  function  of  drag  coefficient:  7  km 
radius  fire;  initial  humidity  ratio  of  7.5  g/kg. 


the  plume  reducing  the  drag.  Thus,  a  realistic  value  for  is  within 
the  range  described.  Even  for  the  lowest  drag  coefficient  value,  the 
displacement  over  12  km  is  more  than  200  m,  i.e.,  more  than  twice  the 
predicted  upper  limit  of  Coriolis  driven  motion  (even  neglecting  the 
countering  effect  of  centrifugal  forces  produced  by  this  additional 
curvature) .  Neglecting  the  Coriolis  force  is  thus  justified  a 
posteriori.  The  effects  of  fire  radius  are  shown  in  Fig.  6  where  the 
results  collapsed  onto  a  single  curve  when  plotted  against  the  product 
of  the  centerline  displacement  (r^,)  and  the  fire  radius  Rf.  This 
allows  extrapolation  to  other  fire  sizes. 


Figure  6.  Plume  centerline  path  for  model  wind  field,  w  =  7.5  g/kg 
and  =  0.5 
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SECTION  3 


THE  EQUILIBRIUM  SMOKE  LAYER 

Next,  we  examine  zone  6  (see  Fig.  1)  and  estimate  the  far  field 
cloud  layer  thickness.  The  estimate  is  independent  of  the  momentum, 
but  depends  on  the  mass  and  energy  in  the  turned  plume.  The  upward 
moving  plume  starts  to  turn  at  an  altitude  H*  where  the  plume  density 
(now  constant)  exactly  equals  the  density  of  the  undisturbed  atmos¬ 
phere  . 

.At  this  altitude,  the  plume  is  neutrally  buoyant  and  has  upward 
momentum  only.  As  a  result,  a  height  overshoot  or  capping  cloud 
develops  (see  Fig.  1).  This  is  a  local  effect;  the  superelevated  part 
of  the  plume  has  a  radius  of  the  order  of  the  plume  radius.  We  con¬ 
centrate  on  determining  the  thickness  and  spread  of  a  flat  layer 
beyond  the  turning  portion  (zone  6),  first  in  a  quiescent  atmosphere, 
and  then  accounting  for  crosswinds. 

3,1  SMOKE  LAYER  THICKNESS  WITH  NO  CROSSWIND. 

The  overturning  process  in  zone  5  mixes  the  plume  gases  resulting 
in  "top  hat"  profiles  of  the  vertical  velocity  and  the  state 
variables.  The  mass  flux  in  the  rising  part  of  the  plume  is: 

M  =  W  ,  (11) 

o  t 

where  P5  =  P4  is  the  plume  density,  Rq  =  ^Rf  is  the  plume  radius  and 
W  the  speed.  This,  in  turn,  is  equal  to  the  flow  at  the  perimeter  of 
the  layer ,  i  .  e  .  . 


M  =  p^  2jrRTU 
o 


(12) 


where  T  is  the  thickness  (at  the  perimeter  of  the  layer)  ,  U  is  the 
speed  in  the  radial  direction,  and  R  is  the  radius.  From  Eqs .  (11) 
and  ( 12) , 
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T 


R^W 

o 


(:l3) 


where  for  steady  flow,  both  U  and  T  are  functions  of  the  radial  dis¬ 
tance  r  only.  The  energy  input  per  unit  time  is 


E 


(14) 


as  the  potential  energy  is  now  zero  relative  to  the  surroundings  (no 
density  difference).  The  rate  of  energy  transported  at  cross  section 
R  is 


+  Mg 


T 

2 


(15) 


The  second  term  represents  the  energy  required  to  maintain  the  layer 
of  smoke  at  a  gi'ten  height.  This  energy  is  calculated  by  ass’iming 
that  the  laver  centered  at  height  Hq  is  compressed  by  the  plume 
outflow.  For  a  Linearly  decreasing  density  atmosphere,  the  change  in 
•••olume  of  the  lower  column  is.  per  unit  radial  distance 


(  lb ) 


The  pressure  applied  per  unit  radial  distance  is  obtained  by  dividing 
the  plume  weight  produced  per  unit  time  by  the  area,  i.e.. 


o 

^  a 
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rR^ 


(17) 


The  work  used  for  compression  of  the  lower  colujnn  is,  per  unit  time 
from  Eqs.  (16)  and  (17) 
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(18) 


P,  A.  -  fig  I 

Equating  energies  and  assuming  negligible  viscous  interactions  (dis¬ 
sipation)  we  have 


2  2 
U  +  g  T  =  W 

Substituting  Eq .  (13)  into  Eq .  (19), 


(19) 


-  w'n  +  g  ^  =  0  .  (20) 

One  can  obtain  U(R)  from  Eq .  (20)  and  then,  using  Eq .  (13)  solve  for 
T(R)  i.e.,  the  layer  speed  and  thus  the  thickness  are  functions  of 
height  and  distance.  Once  T(R)  is  determined,  one  can  generalize  this 
solution  to  take  crosswinds  into  account,  this  being  the  equivalent  of 
having  a  source  in  an  oncoming  flow. 

Equation  (19)  is  a  third  i.iegree  equation  in  U  which  has  one  of 
three  real  solutions:  1)  a  negative  radial  velocity  asymptotically 
approaching  the  value  of  (-W)  as  R  ->•*>;  2)  a  pair  of  solutions  that 
become  real  (and  thus  of  practical  interest)  only  if  R  >  R^-j-  where  R^^ 
is  larger  than  the  original  stem  radius  Rq  (obtained  at  Zq  the  lower 
edge  of  zone  4)  .  At  R  =  R^j-  these  two  solution  branches  have  a  common 
value  of  U,  one  growing  asymptotically  to  the  finite  value  U  -  Wq  as 
R  ->  33 ,  the  other  tending  to  vanishingly  small  velocities  for  R  -*■  =0 
('see  Fig.  7)  . 

The  only  physically  reasonable  solution  is  the  third  one  with  an 
asymptotically  decreasing  speed  and  a  finite  layer  thickness  at  a 
large  distance.  The  fact  that  this  solution  only  occurs  when  R  >  R^-j- 
is  compatible  with  the  definition  of  zone  6  annular  (i.e.,  the  domain 
R  <  R(,j-  is  not  included)  .  The  first  solution  requires  inflow  that  is 
not  consistent  with  a  steady  state  flow.  The  second  solution  results 
in  a  vanishing  thin  layer  in  which  viscous  and  mixing  effects  (which 
were  neglected)  would  become  dominant.  The  third  solution  results  in 
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Figure  7.  Solution  branches  for  outflow  layer  velocity.  Solid  curve  is  the 
only  physically  valid  solution  (Rf  =  7  km,  =  70  m/s). 

a  relatively  thick  layer  in  which  edge  effects  are  small.  An 
asymptotic  cloud  thickness  (see  Fig.  8)  of  around  500  ra  is  obtained 
for  cloud  radii  of  100  km.  The  third  physically  significant  solution 
has  a  domain  of  existence  R  >  Rcr-  The  critical  radius  which  defines 
Che  inner  limit  of  the  horizontal  spreading  zone  is  calculated  as  a 
function  of  the  initial  radius  Rq  =  ^  Rf  and  vertical  velocity  W. 

Figure  9  hows  the  value  of  the  quantity  R^-^/R^ .  The  value  of  the 
radial  velocity  U  at  Che  critical  radius  also  appears  in  Fig.  9.  Both 
velocities  are  linear  functions  of  the  parameter  R(,^/R^  for  a  wide 
inge  of  upward  velocities. 

3.2  LAYER  THICKNESS  WITH  CROSSWIND. 

Numerical  solutions  for  large  fire  plumes  [Bacon  et  al . ,  1987; 

Penner  et  al.,  1986;  Small  and  Heikf>s,  1989]  show  rapidly  decelerating 
vertical  motions  as  the  buoyancy  decreases  and  eventually  turns 
negative.  Thus,  even  in  a  constant  ambient  wind,  the  plume  tilts  or 
shears  more  at  the  higher  altitudes.  This  has  been  observed,  for 
example,  in  Che  Canadian  experimental  fires  [Heikes  and  Small, 
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Figure  8.  Solution  branches  for  outflow  layer  thickness.  Solid  curve  is 
the  physically  meaningful  solution  (Rf  =  7  km,  =  70  m/s). 


Figure  9.  Outflow  velocity  at  critical  radius  as  a  function  of  plume 
centerline  velocity. 


17 


1990b]  .  The  curving  of  the  plume  centerline  is  accentuated  by  the 
fact  that  windspeeds  usually  increase  with  altitude  (see  Fig.  4) .  The 
rising  part  is  in  principle  described  by  circular  cross  sections  in  an 
obliquely  curved  plume.  The  local  angle  to  the  vertical  is  the  ratio 
of  average  vertical  velocity  to  horizontal  velocity. 

An  additional  process  takes  place  at  the  top  of  the  cloud  layer. 
Since  this  layer  is  produced  by  an  overshoot  of  the  plume  gases,  the 
gases  reach  an  altitude  where  their  vertical  velocity  is  vanishingly 
small  before  falling  back  and  spreading  out  radially.  With  no  exter¬ 
nal  wind,  this  occurs  on  the  outside  of  the  ascending  plume.  Here, 
however,  the  windward  part  is  blown  Cowards  the  plume  center  and  is 
shifted  sideways  and  downstream. 

When  the  smoke  layer  reaches  its  equilibrium  height,  the  spread 
can  be  approximated  by  a  simple  two  dimensional  radial  source  flow. 

The  flow  is  deflected  by  the  wind  into  an  elongated  shape,  calculated 
simply  by  superimposing  the  stream  function  representing  the  ambient 
wind  which  is  [Milne-Thorapson,  1968] 

0  =  V  2  (21) 

w  w 


on  a  source  stream  function 


0  =  U  R 

s  cr  cr 


(22) 


where  B  is  the  circumferential  coordinate.  Since  viscous  forces  are 
neglected,  the  combined  stream  function  is 


0  =  0  +0 

W  S 


(23) 


The  cloud  shape  is  thus  a  Rankine  profile  of  finite  length,  obtained 
by  plotting  the  curve  0=0  (see  Fig.  10) .  The  length  of  the  Rankine 
wake  (Xl  in  Fig.  10)  is  a  function  of  the  time  elapsed  since  the  plume 
started  spreading.  By  conservation  of  mass  using  the  layer  thickness, 
the  relation  between  the  time  elapsed  since  plume  inception  (At)  and 
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the  plume  length  is  obtained.  Figure  10  shows  nondimens ional  wake 
length 


as  a  function  of  the  nondimens ional  layer  area  written  as 


(24) 


(25) 


For  a  given  fire,  S'  is  proportional  to  the  time  elapsed  and  the 
length  of  the  layer  grows  almost  linearly  with  time,  as  expected. 


Figure  10.  Nondimensional  length  of  outflow  layer  as  a  function  of 
tne  nondimensional  layer  area. 
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SECTION  4 
CONCLUSIONS 


We  have  shown  that  the  maximum  plume  rise,  altitude  of  the  out¬ 
flow  smoke  layer,  and  cloud  spread  rate  (outflow  velocity)  while 
functions  of  the  prescribed  forcing,  can  be  simply  estimated.  Al¬ 
though  cloud  heights  were  estimated  for  large  area  fires,  a  similar 
analysis  could  be  developed  for  nuclear  dust  clouds. 

The  outflow  layer,  or  mesoscale  smoke  cloud  altitude  and  spread 
rate  are  determined  by  patching  one-dimensional  plume  and  equilibrium 
layer  solutions  through  an  energy  formulation.  The  effects  of 
crosswind  are  then  accounted  for  by  constructing  a  source  function  for 
the  spreading  plume  and  superimposing  the  ambient  winds.  The  solu¬ 
tions,  although  approximate,  provide  rapid  and  reasonably  accurate 
estimates  of  the  smoke  cloud  height  and  the  spreading  cloud 
thickness.  Such  solutions  could  be  used  in  evaluating  mesoscale 
spread  of  smoke  or  dust  layers  and  obscurations  resulting  from  mul¬ 
tiple  nuclear  bursts. 
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